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ABSTRACT 

Context. Weak tensing is one of the best available diagnostic tools to measure the total density profiles of distant clusters of galaxies. 
Unfortunately, it suffers from the well-known mass-sheet degeneracy, so that weak tensing analyses cannot lead to fully reliable 
determinations of the total mass of the clusters. One possible way to set the relevant scale of the density profile would be to make a 
direct measurement of the magnification produced by the clusters as gravitational lenses; in the past this objective has been addressed 
in a number of ways, but with no significant success. 

Aims. In this paper we revisit a suggestion made a few years ago for this general purpose, based on the use of the Fundamental Plane 
as a standard rod for early-type galaxies. Here we move one step further, beyond the simple outline of the idea given earlier, and 
quantify some statistical properties of this innovative diagnostic tool, with the final goal of identifying clear guidelines for a future 
observational test on concrete cases, which turns out to be well within the current instrument capabilities. 

Methods. The study is carried out by discussing the statistical properties of Fundamental Plane measurements for a sample of early- 
type source galaxies behind a massive cluster, for which a weak lensing analysis is assumed to be available. Some general results are 
first obtained analytically and then tested and extended by means of dedicated simulations. 

Results. We proceed to study the best strategy to use Fundamental Plane measurements to determine the mass scale of a given cluster 
and find that the optimal choice is that of a sample of early-type galaxies behind the cluster distributed approximately uniformly in 
the sky. We discuss the role of the redshift distribution of the source galaxies, in relation to the redshift of the lensing cluster and to 
the limitations of Fundamental Plane measurements. Simple simulations are carried out for clusters with intrinsic properties similar 
to those of the Coma cluster. We also show that, within a realistic cosmological scenario, substructures do not contribute much to the 
magnification signal that we are looking for, but only add a modest amount of scatter. 

Conclusions. We find that for a massive cluster (M200 > 10'^ Mq) located at redshift 0.3 ± 0.1, a set of about 20 Fundamental Plane 
measurements, combined with a good weak lensing analysis, should be able to lead to a mass determination with a precision of 20 % 
or better. 

Key words. Gravitational lensing: weak - Galaxies: clusters: general - Galaxies: fundamental parameters - Cosmology: dark matter 

1 . Introduction plane where k{0) + 1 . One can assume that the surface mass den- 

^ , , , , ,. ., . sity vanishes at the boundaries of the image, far from the lens. 

Weak lensmg is a powerful tool to probe the mass distribution j ^^^^ ^^1^^ ^ ^j^^^ ^j^^ ^^^^^ 

of massive clusters of galaxies. Based on the study of the dis- ^^j^ observation is zero. However, this assumption requires 

tortion induced by the lens on images of extended background ^^j^ ^^^^ sufficiently large, which is not always 

sources, weak lensing techniques have been often used to mea- j^j^ Moreover, current structure formation models predict 

d '^''f^l f '^l.';'stersofgalaxies (see e.g. Lombardi et aLM ^ j ^ j ^ ^ nonvanishing surface mass 

■ - - Oowe & Schneiderm Broadhurst et al 20Q1 Clowe et al. densities far from the lens center, so that such an assumption 

m Gavazzi et al.EOOi. However, weak lensmg suffers from ^^^^^ j^^^ ^^^^j ^^^^ estimates significantly underes- 

a fundamental limitation set by the so-called mass-sheet degen- ^.^^^^^ ^^^^^j^^^ possibility is to set = at its minimum so 

eracy . the projected surface mass density map m can be deter- ^^^^ ^^^^ ^^^^ . ^ everywhere positive, which might seem 

mined only up to transformations of the form ^ pl^^^^^j^ assumption. The problem with this approach is that 

;ic(0) i_> a:'(0) = Xk{&) h- 1 - /I. (1) noise can produce negative values of and adjusting the overall 

density profile on the basis of a noise feature may not be wise. A 

In other words, with weak lensing measurements alone it is not popular solution is to fit weak lensing measurements to pa ramet- 

possible to constrain the total mass of a lens, unless further as- ric model mass distributions, usually NFW (Navarro et al. [T997b 

sumptions are made. profiles. This method has the disadvantage that it relies on an as- 

Although many strategies have been proposed to break this sumed mass profile for the lens, which might differ from that of 

degeneracy, no definitive solution has been found so far In prin- the actual cluster, thus reducing the power of weak lensing as a 

ciple, the mass-sheet degeneracy can be removed with the deter- mass measurement technique, 
mination of the absolute value of /c at a single point in the lens 



o 

(N 
0^ 



O 

U 

o 

» 



(N 
> 
(N 

\6 
o 



X 



1 



A. Somienfeld et al.: Direct measurements of the lensing magnification by galaxy clusters 



If we wish to find an assumption-free method to break the 
mass-sheet degeneracy, additional information must be added 
to the weak lensing measurements. Several attempts have been 
made in this direction. One possibility is to include magnifica- 
tion information in the data set. In fact, measurements of the 
magnification and of the shear field can lead to a direct deter- 
mination of the surface mass density a:. Broadhurst et al. (11995) 
proposed a method based on the study of the number counts of 
faint background galaxies for the determination of the magnifi- 
cation. Their technique was successfully used in a few cases of 
particularly massive clusters (Fort et al. 119971 Taylor et al. |1998l 
Broadhurst et al. 2005, Umetsu et al. 20101. However, for this 
method a detailed calibration of nontrivial model quantities, such 
as the number counts of unlensed sources, is essential. Moreover, 
the count process in the central regions of rich clusters is made 
difficult by bright cluster members that hide the faintest back- 
ground galaxies. 

The form of the invariance transformation dl} is referred to 
a fixed source redshift: the same transformation referred to a 
source at a different redshift changes through a coefficient that 
multiplies the term I - A. In other words, each portion of the 
redshift space suffers from a different invariance transforma- 
tion, so that in principle the mass-sheet degeneracy can also 
be broken by combining lensing information from images of 
sources at different known redshifts. Bradac et al. (l2004i inves- 
tigated this possibility in the context of weak lensing, consider- 
ing the hypothetical case in which the individual redshifts of the 
lensed background galaxies are available. They showed that, un- 
der these favourable circumstances, the mass-sheet degeneracy 
can be broken for critical clusters (i.e. those that can produce 
multiple images), but still not for subcritical ones. 

A significant improvement can come from the addition of in- 
formation from strongly lensed images (i.e. arcs or multiple im- 
ages), provided that the redshifts of the strongly lensed sources 
differ from the mean redshift of the sources used for the weak 
lensing analysis. Indeed, several methods based on the inclu- 
sion of strong lensing information have been devised and ap- 
pUed to real cases (Bradac et al. '2005a,' '2005b,' Cacciato et al. 
120061 Diego et al. 2007, Merten et al. 2009). However, strongly 
lensed images can only be produced by critical clusters. Thus, 
for subcritical clusters the mass-sheet degeneracy still remains a 
fundamental problem in the determination of the total mass. 

It is in this context that Bertin & Lombardi (2006, BL06 from 
now on) proposed a new method to measure the lensing magni- 
fication induced by a cluster, which can be used to break the 
mass-sheet degeneracy. They showed that estimates of the mag- 
nification can be obtained by observing background early-type 
galaxies. Early-type galaxies can be treated as standard rods, in 
virtue of the empirical law of the Fundamental Plane. From a 
Fundamental Plane measurement, the intrinsic effective radius 
of an early-type galaxy can be recovered with a 15% accuracy 
( IBL06I ). Then, by measuring the observed (magnified) effective 
radius, the magnification can be derived. 

In this paper we investigate further the possibilities opened 
by this technique. In particular, we wish to determine which ac- 
curacy on the determination of the total mass of a cluster can be 
achieved by combining Fundamental Plane measurements with 
a weak lensing study. This is done by prescribing a method to 
break the degeneracy and by studying its properties, both analyt- 
ically and with the aid of simulations. In principle, substructures 
in the lensing cluster can introduce noise in the magnification 
signal, so that a given measurement might be used to set interest- 
ing constraints on the amount of substructure rather than on the 
mass of the cluster. This issue is also studied in this paper, with 



the use of numerical simulations. Then we address the problem 
of identifying the optimal conditions on lens and source-galaxies 
in order to get a significant measurement in concrete cases. 

The structure of the paper is the following. In Sect. |2] we 
give the basic lensing equations and present the problem of the 
mass-sheet degeneracy. In Sect. [3] we show how Fundamental 
Plane measurements can be used to infer the magnification and 
present a simple method to use such measurements to break the 
mass-sheet degeneracy. In Sect. |4] the statistical properties of 
this method are studied. In Sect. |5]we address the issue of how 
substructures can influence the magnification signal we seek to 
measure. In Sect. |6] we describe the simulations set up to test the 
method and show the results. Conclusions are drawn in Sect. [T] 



2. Weak lensing preliminaries 

2.1. Basic notation and equations 

We start by introducing the projected surface mass density of 
a given lens, 1.(9). The nature of the lensing equations make it 
convenient to introduce the dimensionless surface mass density 
k{0, z) for a source at redshift z, defined as 



K{e,z) = 



Icriz) 



where S^r = 



4nG DdDds 



(2) 



Here Dd, Ds, Dds are angular diameter distances of the lens and 
source with respect to the observer, and of the source with re- 
spect to the lens, respectively. The surface mass density is re- 
ferred to a fiducial source at infinite redshift. 



K{0,z) = z(z)K(e), 

through the cosmological weight function 
Z(z) = — — H{z - Zd)- 

Icriz) 



(3) 



(4) 



Here H(z-Zs) is the Heaviside step function, to take into account 
that images of sources at redshift lower than that of the lens are 
not lensed. 

An important quantity that enters the weak lensing problem 
is the redshift-dependent reduced shear g(0, z): 



8(0, z) 



Z(z)y(0) 
l-Z(z)K(0y 



(5) 



where the shear y{0) can be expressed as a nonlocal function of 
K{0): 



y{0)^- { D{0-0')K{0')d^9', 
1 



with D(0) = - 



(6) 



(7) 



Weak lensing consists in the study of the distortion induced by 
the lens on the images of background sources. This can be done 
in terms of a complex ellipticity e, defined from the quadrupole 
moments Qij of the surface brightness as 



611-622+2/612 



611 + 622+2^611622-6 



(8) 
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Seitz & Schneider d 19971 1 showed that the image elHpticity is 
related to the intrinsic (unlensed) source ellipticity e' in the fol- 
lowing way: 



if\8(0,z)\ < 1 



otherwise 



(9) 



Therefore, under the assumption that the intrinsic ellipticity dis- 
tribution of background sources is isotropic, the expectation 
value of the observed ellipticity for sources at redshift z is: 



E[e(z)] = 



1 



if\g(0,z)\<l 



otherwise 



(10) 



In the general case of sources distributed in redshift, the follow- 
ing approximation holds for k < 0.6: 



E[e] 



{Z}y(0) 



(11) 



as shown by Seitz & Schneider ("1997). (Z") are the moments of 
the redshift probability distribution of the background galaxies. 

2.2. The mass-sheet degeneracy 

By averaging over image ellipticities of background galaxies and 
identifying (e) with E[e], we can estimate the quantity (fTTT i in 
the field of observation. Then, since y depends on k through Eq. 

relation (fTTT i can be inverted and the surface mass density k 
can be recovered from the observed average ellipticity. Practical 
realizations of this picture were provided by Kaiser (119951 . Seitz 
& Schneider (1997) , Lombardi & Bertin (1999 l. 

It can be shown that the quantity E[e] in Eq. (fTTT i. which is 
the observable quantity, is invariant under transformations of the 
form 



k(0) k'{0) = Ak{0) + w{\- A), 
where we introduced 

<Z) 

W = 



(12) 



(13) 



This is the mass-sheet degeneracy for the general case of sources 
distributed in redshift: with ellipticity measurements it is only 
possible to recover the surface mass density up to the above 
transformation. In the simplified case of sources all at the same 
redshift z, the invariance transformation reduces to ([1}, with 

K(e,z) = K(0). 

In principle, the mass-sheet degeneracy can be broken with 
a local measurement of the magnification. In fact, the magnifi- 
cation jj(0, z) is related to k(0) and y(0) through the following 
relation: 



ti(0,z) = [1 - z{zM0)f - zHzM0r 



(14) 



By measuring (e) and yu(z) and by combining (fTTT i with (fT4T l we 
obtain a two-equation system in terms of k and y, which can 
be solved for k, thus breaking the mass-sheet degeneracy. This 
requires that we know which redshift z the magnification mea- 
surement is referred to, in order to calculate the cosmological 
weight Z(z) that enters Eq. (fT4T l. 



3. Breaking the mass-sheet degeneracy 

In this section we will introduce a mass measurement method 
based on the combined use of weak lensing and magnification 
measurements. 



3.1. The Fundamental Plane 



The Fundamental Plane (FP from now on; Dressier et al. 119871 
Djorgovski & Davis 1987) is an empirical scaling law that ap- 
plies to early-type galaxies (E/SO). It relates three well-defined 
observable quantities for these objects: the effective radius, Rg, 
the effective surface brightness, (SB)^, and the central velocity 
dispersion of the stellar component, ctq. The three quantities are 
related in the following way: 

Log/;, = Log r, + LogDA(z) = o-Logcro H-/3<SB), + y, (15) 

where a, p, and y are empirically determined coefficients that de- 
pend on the waveband of observation, rg is the effective radius in 
angular units and Da{z) is the angular diameter distance of the 
galaxy at redshift z. The existence of such a relation has been 
extensively confirmed by a number of studies on both field and 
cluster galaxies out to cosmological distances (e.g. J0rgensen et 
al. fT993T Bender et al. fT998l Treu et al. 1999 1. The measurement 
of the Fundamental Plane parameters for the most distant sam- 
ple of objects has been carried out by van der Wei et al. (120051 
vdW05 from now on), who examined early-type galaxies out to 
z ~ 1.1. This relation is observed to hold within a 0.07 scatter 
on Log re, or 15% on r^, rather independently of the position on 
the FP plane (J0rgensen et al. il996) and increasing with increas- 
ing redshift (Treu et al. I2005I I. Treu et al. (12005) have quanti- 
fied the increase in the scatter in Log rg of the FP relation as 
dcTy/dz = 0.032 + 0.012, which translates into a scatter of 23% 
in at z = 1. It is still not clear if the source of this scatter 
is totally intrinsic or if it can be reduced by improving the ob- 
servational precision. Auger et al. (120 lOt estimated the intrinsic 
scatter of the FP to be as low as 11%. 

Observations have pointed out a variation of the coefficient 
7 with redshift, quantified as dy/dz = 0.58^QQg by Treu et al. 
(I2005I I. For cluster galaxies there is evidence for a slower evolu- 
tion (see e.g. Wuyts et al. 2004 vdW05), so that the Fundamental 
Plane for cluster galaxies appears to differ from that of field 
galaxies. In vdW05 it was shown that this difference is not signif- 
icant for massive (M > 2 x lO^M©) galaxies, and a similar result 
was found by van Dokkum & van der Marel (I2007t . According 
to vdW05, the scatter of the FP is also smaller for the more mas- 
sive objects. Evidence for a variation of the coefficients a and 
with redshift has also been reported (Treu et al. I2005K but this is 
generally taken to be less significant. 



3.2. The FP seen through a lens 

As shown in BL06, the Fundamental Plane changes in a well 
defined way when viewed through a gravitational lens. Both the 
surface brightness (SB), and the central velocity dispersion ctq 
are lens-invariant. Therefore, by measuring these two quantities 
for early-type galaxies and by making use of the Fundamental 
Plane relation (fTsT i it is possible to obtain an estimate of the ef- 
fective radius of the observed galaxy, Rf^\ By measuring the 
redshift of the galaxy it is then possible to convert this measure- 
ment of the effective radius to angular units, rf^K By observing 
the effective radius 4°*"*', magnified by the lens effect, the lens 
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magnification will then be given by the square of the ratio of the 
intrinsic size to the observed image size, as 



,(FP) 



(16) 



With a 15% scatter of the Fundamental Plane relation, the same 
error will affect our estimate of the intrinsic effective radius, 
rf^^ , which translates into a ~ 30% error on the magnification. 

Actually, the quantity that is lens invariant is the intrinsic sur- 
face brightness, which in general differs from the observed sur- 
face brightness. Nevertheless, van der Wei et al. (2005 1 showed 
that the intrinsic surface brightness of galaxies at high redshift 
can be effectively measured by fitting Sersic models convolved 
with the PSF. We expect this task to be made easier by the mag- 
nifying effect of lensing. 

3.3. A minimum-x^ approach 

The method presented here has been developed and tested for 
noncritical lenses only, although it can be generalized to the crit- 
ical case. Therefore, from now on we will assume that the lens 
does not have critical curves, unless stated differently. This is 
also the most interesting case, because it is for subcritical lenses 
that the problem of the mass-sheet degeneracy is harder to over- 
come. 

Suppose we performed a weak lensing analysis of a cluster, 
which led to the determination of the average distortion (e) {6) 
within the field of view. By inverting the distortion map (for 
instance with the nonparametric method of Lombardi & Bertin 
11999 ) it is possible to recover the surface mass density of the lens 
k(6) up to the invariance transformation (fTZt . Then, suppose that 
we performed a set of A^fp Fundamental Plane measurements, 
through which we estimate the magnification in a corresponding 
number of positions {0, ) on the image plane. 

Among the infinite possible mass density maps compatible 
with the distortion measurements, spanned by (fTZt . we select 
the one for which the accordance with the observations of mag- 
nied early-type galaxies is best. To do so, we first transform 
the estimates of the magnification into estimates of the surface 
mass density, K^^^\9i), making use of Eqs. (fTTl l and (fT4l l. This re- 
quires that we know the average distortion in the image position, 
(e) (Oi), determined from the weak lensing study, and the redshift 
Zi of each galaxy, which can be measured with the same spectro- 
scopic observation necessary to measure ctq. The expression that 
gives K^^\6i) in terms of /i(0,), (e) (0,) and Zi is the following: 



.(FP) 



-b- V^2^-fl(c^^I/^ 



where 



wHzy 



\{e}\ 



b = Zi 



w{zY 



c= 1 



{zy 



and Z; = Z(z,). 



(17) 



(18) 



(19) 



The minus sign for the square root in (fTTt comes from the fact 
that the lens was assumed to be everywhere subcritical. 

The fit is performed by minimizing the following penalty 
function: 



X 



(20) 



where kq{6) is the surface mass density distribution inferred from 
the weak lensing reconstruction, and cP: are suitably chosen 
weights. The minimum-;^^- condition can be found by imposing 
dx^/dA = 0, which leads to the determination of the estimator A 
for the parameter A: 



y ^[^r'-w][^o(e,)-w] 

A'fp 1 

2— [^o(0')->^]' 



(21) 



The choice of the weights erf that enter this function requires 
particular care. Generally, in minimum-;t'^ approaches the values 
of cr, are taken to be proportional to the measurement errors of 
the quantity over which the fit is performed: in the present case, 
the surface mass density k^^^\ By assuming that the errors come 
only from the Fundamental Plane measurements (i.e. under the 
assumption of perfect weak lensing measurements), error prop- 



agation on K^^^^ gives 



[l-Z(z)Kf-^{l-^X\{e)\H0) 
(ZY \ wl 



A/Cfp = crFp(z)- 



(22) 



Z(z)[l - Z{z)k\ - (l - -] I {s) \\e) 

where crFp(z) is the scatter in of the Fundamental Plane. 
Therefore, A/cfp depends on the value of k, which is the quan- 
tity we are trying to determine with the fit. This complicates 
the definition of the weights. One would be tempted to define 
CT; = AKi/^HOi)), but by doing so we would introduce a bias in 
the estimate of A: this is because a higher weight would be given 
to measurements in which fluctuations give higher values of the 
magnification, relative to cases in which ^ is underestimated. 
Then the fit would yield preferentially higher mass estimates. 

The solution we propose is an iteration procedure, seeded by 
the definition of the weights based on the model surface mass 
density kq obtained from weak lensing: in (l22T i we identify kq 
with K and define 



CTi = AKpp{Ko{9i)). 



(23) 



We can then minimize the;i^^ function and obtain a new estimate 
of the density map. If Aio) is the value of the parameter A obtained 
from the fit with ( 1211 1. then this new model for the surface mass 
density will be given by 



Kl = A(0)Ko + W(l- /l(())). 



(24) 



At this point we can update the weights to this new model by 
redefining cr, = AKpp[Ki(Oi)] and iterate the procedure. Fig. [1] 
outlines the basic steps of the method. The method converges 
quickly and the final solution is found to be invariant under trans- 
formations of the form (fT2t of the input surface mass density 
map K(). 

Eq. ( l20l i is only one possible definition of a x^ function for 
our problem. We could have defined a x^ in terms of other quan- 
tities, such as the magnification /i or its inverse | detA|. The ad- 
vantage of our choice lies in the fact that the minimum condi- 
tion oljl^V^'l = is a linear function of A. This makes it relatively 
easy to study the statistical properties of the method, which will 
be the subject of Sect. |4] 
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Weak lensine reconstruction 



where we introduced 



n = 



Iteration 
n = n + \ 



a-,=^KfAK'"\\{B)\,Zi) 



Fit 



Fig. 1. The fitting process. 



4. Statistical properties of the method 

To assess the reliability of the method introduced to break the 
mass-sheet degeneracy, it is necessary to study in detail the sta- 
tistical properties of the results obtained. In particular, we wish 
to clarify if the estimate of the total mass resulting from this pro- 
cedure is biased, and to determine the expected uncertainty. 

In order to do this, we need to discuss the conditions of ob- 
servation first. Obviously, it would be preferable to work with 
a large number of measurements, but a practical limit is im- 
posed by the telescope time required to measure velocity dis- 
persions. We believe that a reasonable number of Fundamental 
Plane measurements that can be efifectively carried out in an ob- 
servational campaign is A^pp a: 20. Another important issue is 
the possibility of nding a sufficient number of sources behind the 
observed gravitational lens. BL06 estimated the number density 
of practically observable sources to be ~ 2arcmin"^. This es- 
timate was based on a study of magnitude-number counts of 
early-type galaxies (Glazebrook et al. 1995 ). This sets a limit on 
the field of view required to have a sufficiently high statistics of 
Fundamental Plane measurements. It also tells us that, unless a 
large field of view is examined, the observer does not have much 
freedom in selecting objects. Our statistical study will take into 
acocunt these important factors. 

In this section, the weak lensing study will be assumed to 
provide a perfect reconstruction of the surface mass density map, 
except for the mass-sheet degeneracy. Since the fit method is in- 
variant under tranformations of the form (ITZl on the input den- 
sity map kq{6), we can also assume without further restrictions 
that kq{0) is the exact surface mass density of the lens: 



(25) 



Based on these assumptions, we will study the statistical proper- 
ties of the estimator A at the first iteration (i.e. /I(0)). The results 
obtained will provide information on the error in the determina- 
tion of the mass density map KiiO) after the first iteration. Since 
the method converges quickly, there is little difference between 
K\ (0) and the final density map, therefore we will consider the 
difference between ki and the exact density map k to evaluate 



— [^r^-w]we,)-w] 



(27) 



We assume that the number of Fundamental Plane measurements 
Npp is fixed, and that the image positions are randomly dis- 
tributed. We will denote by Pe~{9, z) the probability distribution 
in redshift and image position of the observed galaxies. 

Given these definitions, we write the quantities E{/,) that en- 
ter (l26l l as 

E{/,) = ]~[ I PgAOj,zj)d0jdzj I PMr^)dKf^hi, (28) 

where is the probabiUty distribution for the estimates of k 
from a Fundamental Plane measurement. Let us make the as- 
sumption that Pit{Kf^^) is a Gaussian function with dispersion 
cr,, centered on the exact value of the surface mass density: 



P.[Kf\{K,a,]) 



1 



exp 



„(FP) 



K{0>)V 



(29) 



The validity of this assumption will be discussed in Appendix A. 
We then perform the integration in dK^f^^ in (l28T l, which yields 



dK] 'P^(k) 



— HOd-wY 



(30) 



from which, together with 

Eai-i. 



it easily follows that 



(31) 



Since the input density map was assumed to be the exact surface 
mass density of the lens (/fo(0) = k(0)), this result implies that 
the expectation value for the reconstructed surface mass density 
distribution is the exact one: the estimator A is not biased. 

As a second step, we study the second moment of the proba- 
bility distribution for A, the variance: Var(/l) = E{A-} - E^{A]. A 
similar calculation leads to 



Var(i) = E{A^} - E^{A} = 



1 



—[K(0d-wr 



(32) 



where angle brackets indicate that the expression must be aver- 
aged over the possible image positions and source redshifts. 

We note that in the limit of high A^fp the variance has the 
typical behavior ~ 1 /Npp. We have exactly Var(l) oc 1 /Npp in 
the case of a uniform sheet of constant surface mass density. 



4.1. Error on A 

We define the expectation value of the estimator A as 

Ea) = 2E{/,), 

1=1 



(26) 



4.2. Error on the total mass 

We have just shown the basic results on the accuracy of the 
method in determining the parameter A. Now we will examine 
the problem of how an error on A translates into errors on the 
estimated total mass M of the lens. 
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The total mass inside the field of view is given by the integral 
of K over the observed poition of the image plane 0: 



J© 



where Dd is the angular diameter distance of the lens relative to 
the observer The surface mass density is subject to the invari- 
ance transformation (fTZt . The expectation value of the measured 
mass M is 



as a function of /c (i. e. image position) for a source at redshift 
Zs - 0.6, 0.8, 1.0. For this and the following examples, the value 
of w is computed by assuming a redshift distribution of weak 

(33) lensing measurements P'"^^\z) z^exp(-z/zo) with zo = 2/3, 
as in Bradac et al. (I2004I I. while the scatter of the Fundamental 
Plane is assumed to be ctrp = 0.15. It can be seen that Ji in- 



E{M) 



Ak{0) + w{\ - A)\(t0 \ ^ M, 



(34) 



where we made use of dSTJ, and Zee = ^criZs — > °°)- If the es- 
timator for A is not biased, then there is also no bias on the 
estimate of the total mass. 

Let us consider the variance of M. A straightforward calcu- 
lation based on Eq. ( 1311 1 gives 



Var(M) = E{M'-] - E^jM) = M - wE, 



rDl f Si 
Je 



0\ Var(i).(35) 



By taking the square root of (135) and dividing it by M we obtain 
the expected relative dispersion of M: 



o-(M) 
M 



1 



(T{A), 



(36) 



where k is the average surface mass density inside the field of 
view and a-{A) is the square root of (l32l i. namely the dispersion 
of 1 

This result allows us to calculate, for a given lens, the ac- 
curacy with which its mass can be measured. For fixed average 
surface mass density k, the only quantity that determines this 
accuracy is the precision in the determination of A, a-(A)- This 
quantity might depend on the shape of the lens mass distribu- 
tion, so that lenses with certain characteristics may be more suit- 
able candidates than others. This issue, which is important for 
determining which are the ideal lens candidates for an applica- 
tion of the present technicjue, is addressed in the next subsection 
by studying o-{A) and cr(M) for different lens models. 

4.3. Simple examples 

To better understand the above results, it is useful to focus on 
simple situations. We start by discussing the case with A^fp = 1 
and consider the quantity in brackets in (l32t : 



Ji = 



1 



(37) 



Then, as a first example, we consider an approximate represen- 
tation of the Nonsingular Isothermal Sphere (NIS), with density 
profile given by 



1 



The projected mass density K(ff) of a NIS is given by 
1 1 



(38) 



(39) 



We take cr,, = 1000 km s" and = 58.8 kpc and place the 
lens at redshift Zd = 0.3, so that the projected surface mass den- 
sity at the center is k(Q) = 1. In Fig. |2] we plot the quantity Ji 




Fig. 2. Ji, as defined in ( l37l l. as a function of k (i.e. image posi- 



tion) for a NIS lens with cr,, = 1000 km s 
for three source redshifts. 



and Tc = 58.8 kpc. 



creases with increasing k, although mildly. This result suggests 
that, for a given lens observed within a given field of view, it 
would be preferable to perform Fundamental Plane measure- 
ments on objects whose images lie where the surface mass den- 
sity is lower, as the expected dispersion on A is lower. However, 
it is important to recall the assumptions that underlie this result. 
In particular, we are assuming that no error comes from the we ak 
lensing analysis. When this assumption is dropped the situation 
changes. Real cases are more complex: as will be shown in Sect. 
6, weak lensing reconstructions tend to underestimate the sur- 
face mass density in the central parts of lenses and overestimate 
it in the outskirts. Therefore, normalizing the overall mass scale 
of a lens with Fundamental Plane measurements limited to par- 
ticular regions of the image plane can be risky, because biases 
can be introduced. Weak lensing errors can be better tackled by 
sampling the image plane uniformly. This issue is discussed fur- 
ther in Sect. 6.2. 

As a second step, we will consider the variation of Var(/i) 
as a function of the average surface mass density k of the lens, 
for two simple lens models, again in the case A^fp = 1 . As lens 
models we consider NIS lenses with various values of cr,, and 
but fixed central surface mass density /Cc = 1 . Equation ( l39l l 
shows that NIS lenses with different values of r^ and fixed k^- 
are rescaled just in the angular dimension. Therefore, for a given 
lens, changing r^ and keeping Kc fixed is equivalent to choosing 
a different field of view. 

For comparison, we also considered lenses of constant sur- 
face mass density K(ff) - k and no shear In Fig. [3] we plot the 
variance of /I as a function of k for these two lens models, for 
a single source at fixed redshift Zs - 0.8. The space average in 
Eq. ( l32l l was calculated by ignoring effects of the magnification 
on the image position probability PgiO). The value of the vari- 
ance of the estimator A is found to be similar for the two kinds 
of lens, for fixed average density k. This result is presumably a 
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K 

Fig. 3. Var(A) as a function of the average surface mass density k 
for two lens models: a NIS with /c(0) - 1 (solid line) and a sheet 
of constant surface mass density and no shear (dotted line). 



consequence of the mild dependence of Ji on k, as observed in 
the plot of Fig.|2] 

Finally, we calculate the dispersion on the measurement of 
the mass, o-{M)/M for the two cases considered above, plotting 
the results in Fig.|4] The results for the two lenses are practically 
indistinguishable. This result, which reflects the similarity of the 
values of o-(A) found for the two lenses, suggests that the shape 
of the mass distribution plays little role in determining the accu- 
racy of the mass measurement, while the decisive factor is the 
average surface mass density within the field of view, k. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0. 



K 



Fig. 4. o-(M)/M as a function of the average surface mass density 
K for two lens models: a NIS with k{Q) - 1 (solid line) and a sheet 
of constant surface mass density and no shear (dotted fine). 



already mild dependence of Var(/i) on the image position (see 
Fig. IS, we do not expect our results to change significantly. 



5. Substructure effects 

In some galaxy-galaxy strong lensing systems it has been noted 
that while it is relatively easy to build smooth lens models that 
reproduce well the multiple image positions of distant QSOs, the 
same models are unable to reproduce the observed flux ratios of 
these images (e.g. Kent & Falco [T988l Hogg & Blandford lT994l 
Falco et al. 1 1 9971 ). It is now common belief that these flux ra- 
tio anomalies may be due to the presence of substructures, and 
several attempts have been made at providing suitable lens mod- 
els that take substructure into account, such as those by Mao 
& Sch neider (1998), Metcalf & Madau ( l200TT l. Bradac et al. 
(120021) . Chiba (120021 ). Their argument is supported by the fol- 
lowing reasoning: a substructure having an Einstein radius com- 
parable to the angular size of the source's light emitting region 
can cause a significant change in its apparent size, while its posi- 
tion on the lens plane can be little affected. This condition is rel- 
atively easy to obtain for the case of a compact source such as a 
distant QSO being lensed by a galaxy; in fact, objects as massive 
as a typical globular cluster (~ 10'' Mq) can change appreciably 
the flux received from such a source, if properly aligned. This 
effect is enhanced in the proximity of critical curves (i.e. curves 
where | detA| ^ 0). Since 



fi = IdetAf' = 



1 



(40) 



if I detA| ^ a small change in k can produce a large change in 
the magnification. 

Weak lensing is a good tracer of the surface mass density 
averaged over finite portions of the image plane, but is not sensi- 
tive to small scale variations of the projected mass distribution. 
In fact, weak lensing methods recover the reduced shear g by av- 
eraging the distortion signal over a number of background galax- 
ies over angular scales of tens of arcseconds (e.g. see Lombardi 
et al. 120001 ). Therefore, they only provide smoothed mass den- 
sity profiles. If we want to break the mass-sheet degeneracy 
with magnification measurements, we must be sure that these 
magnification measurements also reflect the properties of this 
smoothed mass profile. Substructures modify the lensing signal 
in a nontrivial way and can complicate the interpretation of mag- 
nification measurements. For the purpose of constraining the to- 
tal mass of the lens, the contribution from smaller clumps has 
the same eff'ect of noise. 

The problem is addressed with the aid of numerical simu- 
lations. We construct a model cluster as a superposition of a 
smooth principal halo and a number of smaller subhalos. We 
generate images of background early-type galaxies, compare the 
observed magnification with that inferred from a smooth lens 
model, and then analyze the results. 

In this section we will use the terms clump, halo, subclump, 
subhalo, substructure synonymously to refer to mass concentra- 
tions inside a galaxy cluster 



We recall that these results were obtained by assuming that 
the image positions are distributed randomly. This may not be 
true if the observer plays an active role in the object selection. 
For example one may wish to sample the field uniformly, avoid- 
ing close pairs of images. In this way the image positions will be 
correlated. This will aff'ect the value of o-(A). However, given the 



5.1. Modelling substructure in clusters of galaxies 

It is not clear how much substructure is present in clusters of 
galaxies. Cosmological simulations predict the existence of a 
large number of subhalos of mass M < IO'^'Mq, but observations 
have failed to prove their presence, so far. Here we will adopt a 
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conservative approach and will take into account the possibil- 
ity that substructures are present in the abundances predicted by 
ACDM models. 



5.1 .1 . Masses and spatial distribution of subclumps 

N-body simulations of structure formation at cluster scales based 
on the ACDM scenario have shown that the dark matter halos 
that maintain their identity after the formation process can be 
approximately described by a power-law mass function: 



dN 
MI 



oc M 



(41) 



with a ~ 0.9 H- 1, as demonstrated for example in the papers by 
Tormen et al. (1998 ), Ghigna et al. (l2000l l. De Lucia et al. (l2004l l 
and Gao et al. (120041 1. These dark matter subclumps are believed 
to account for a few percent of the total mass of the cluster Note 
that the mass fraction in substructure also depends on the age of 
the cluster: as a cluster evolves the subhalos merge to the main 
halo until the whole structure virializes and fewer substructures 
are left. 

Given these results, we will model clusters in a semi-analytic 
way. A smooth mass density profile for the main halo is chosen; a 
distribution of subclumps is randomly generated from the mass 
distribution given by (l4Tl i. with the total mass in substructures 
fixed; these subclumps are randomly distributed with a spatial 
probability distribution proportional to the mass density of the 
main halo. More details about the practical realization of this 
procedure are given below. 

Numerical simulations have also shown that more massive 
clumps tend to be located far from the cluster center, where 
only small scale halos survive the merging process (Tormen et 
al. 119981 Ghigna et al. I2000I I. For simplicity we adopted a uni- 
form spatial distribution of subhalos, as appropriate for the kind 
of study we wish to perform. 

5.1.2. Internal structure of subclumps 

For the description of the internal structure of dark matter subha- 
los we follow the work of Metcalf & Madau ((2001)), who devel- 
oped numerical simulations to study the effects of subclumps in 
galaxies on the measured lensing magnification of distant QSOs. 
In our work we extend the use of their tools to galaxy clusters 
environments. Metcalf & Madau modeled subclumps as trun- 
cated singular isothermal spheres (TSIS). The advantage of us- 
ing singular isothermal spheres is that their lensing properties 
can be easily described analytically. On the other hand a singular 
isothermal sphere has infinite mass. For this reason a truncation 
radius is introduced. For a clump of mass m at radial distance R 
from the center of the main halo the truncation radius is taken to 
be equal to the tidal radius r,, which is estimated as (Metcalf & 
Madau l200lT) 



3M(R) 



1/3 



(42) 



where M{R) is the mass of the main halo enclosed by the sphere 
of radius R. With this choice, if M(R) grows less steeply than 
R^ (both SIS and NFW models have this property) then clumps 
closer to the center tend to be more compact than those that lie 
far from the center 

Once the mass and truncation radius of the clump are fixed 
there is a unique truncated singular isothermal sphere with those 



characteristics. In particular, if we adopt the following notation 
for the density of a TSIS, 



P(r) 



2nGr^- 




if r < 



elsewhere 



(43) 



the parameter cr,, (often called the velocity dispersion) is given 
in terms of the mass m of the clump by 



= — = — [3M(«)]'/^ . 



2r, 



2R 



(44) 



5.2. Lensing simulations 

5.2.1. General prescriptions 

A simulated cluster is generated with the procedure described in 
Sect. 15.1! For the main halo we adopt a nonsingular isothermal 
sphere (l38T l. Subhalos with a mass distribution given by (l4ll and 
cumulative mass that accounts for a fraction / of the total mass 
of the cluster within r20() are then added to this smooth compo- 
nent. Here r2(K) represents the radius of the sphere whose mean 
density is 200 times the critical density at the redshift of the ob- 
ject, and M2()o is the corresponding mass. The total mass in sub- 
halos is then /M2oo- A direct use of ( 1411 1 results in a very large 
number of clumps with very small mass. These small clumps are 
not relevant for the lensing problem because they produce neg- 
ligible effects on the magnification of extended images such as 
those of early-type galaxies. Moreover, the inclusion of a large 
number of clumps increases the computational effort required to 
run the simulations, so that it is useful to introduce a cutoff on 
the lower range of possible masses. A reasonable choice for this 
lower mass cutoff is a value ot„„„ for which the typical lensing 
deflection angle is only a small fraction of the angular size of the 
source considered. 

To test how the results of our simulations depend on the as- 
sumed internal structure of clumps we also adopted an alterna- 
tive (and non realistic) internal profile: we assumed clumps to be 
point masses. Clearly, the TSIS clump will produce a smaller de- 
flection angle up to r, . For distances larger than r, the two clumps 
will produce equal values of the deflection angle, as they both act 
as if they were point masses located at the center For this rea- 
son, the effects on magnification produced by a point mass will 
be generally higher than those produced by an extended lens of 
the same mass. Therefore, if we study the problem of the mag- 
nification induced by substructure by modelling substructure as 
point masses we can obtain an upper limit on the effects of sub- 
structure on the magnification. 

5.2.2. Practical realization 

We simulate the observed image of a circular early-type galaxy 
with effective radius R^ - 5 kpc at redshift Zs = 0.8 being lensed 
by a NIS cluster with substructures at redshift Zd - 0.3. The pa- 
rameters of the NIS model chosen for the test of the method are 
taken from a study of the Coma cluster by De Boni & Bertin 
(I2008I I. who found the following best fit parameters for the de- 
scription of the dark matter halo as a NIS model: = 88 kpc 
and cr,, - 1 156 km s"'. 

We fix the mass fraction in substructure / and generate sub- 
structures with the procedure described above, with a lower mass 
cutoff of m„,j„ - lO^M©. This value was chosen because, for the 
redshift configuration of our system, the Einstein radii of point 
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masses less massive than m,„,„ are smaller than 0.3 kpc, thus with 
little impact on the distortion of the images chosen for our study. 
Note that in a situation in which all the substructures account 
for a fraction / of the total mass of the cluster, those with mass 
greater than m,„,„ will in general account for a lower fraction 
/' of the total mass. A fraction / - /' will consist of clumps 
with m < m,„i„. Since the lensing effect of these low mass mass 
clumps is small we simulate them by adding a smooth compo- 
nent for a fraction / - /' of the total mass. We generated clusters 
with values of / equal to 0.05, 0.10 and 0.15. In these three sit- 
uations and with the adopted value of m,„,„ the corresponding 
value off is 0.029, 0.061 and 0.093 respectively. 

To simulate the image formation we proceed as follows, tak- 
ing inspiration from the work by Metcalf & Madau (1200 H . We 
define a field of view of 4 x 4 arcmin^ centered on the main halo. 
We then define a grid over the entire field of view. The resolu- 
tion of the grid is such that the observed effective radius for the 
early-type galaxy sources in the absence of magnification is 10 
grid cells long. 

At every grid point the lensing deflection angle is computed. 
The contribution of the subclumps is calculated in the following 
way: the clumps are first assumed to be point masses and the 
deflection angle generated by each clump is calculated at every 
grid point. Then, only for the clumps that lie within the field of 
view considered, the deflection angle inside the circle of radius 
r, is corrected with the expression for the deflection angle of a 
TSIS lens (Metcalf & Madau l200TT) : 



a(x) — ao < 



— - 1 -I- arctan \ —z - 1 if a < 1 



if fl > 1 



(45) 



where a = xrdr,, x = rjvc. More simply, a - rjr,. Differently 
from the work of Metcalf & Madau (2001), we calculate the de- 
flection angle once and for all the grid points in the field, taking 
into account the contribution to a of every subhalo in the lens. 

Once the map of the deflection angle over the field of view 
is created, images of early-type galaxies are generated with a 
uniform distribution on the lens plane. 

Effects of the magnification on the image position distri- 
bution are neglected, i. e. a uniform distribution in the image 
plane is adopted. A rigorous treatment of this aspect would re- 
quire knowing the luminosity function of our target galaxies, but 
this is beyond the goals of this paper Nevertheless, we explored 
different scenarios by applying importance sampling to the out- 
come of our simulations. In particular, we examined two differ- 
ent probability distributions: proportional to the magnification, 
and inversely proportional to the magnification. The change in 
the results with respect to the uniform case was less than 15%, 
indicating that our approximation does not alter the conclusions 
of our study. 

Sources are described with circular de Vaucouleurs profiles, 
truncated at r^,. Images are then constructed in the following way. 
We define a square in the image plane centered on the adopted 
location for the image center The size of this square is chosen 
so that the length of its sides are a few times the effective radius 
times the linear stretching (~ /i'^^) produced by the smooth com- 
ponent of the lens. In this way the observed image is guaranteed 
to lie within this square. Then, each pixel inside the square is 
mapped to the source plane with the lens equation of the sim- 
ulated cluster, where the deflection angle was calculated with 
the procedure described above. Since gravitational lensing pre- 
serves surface brightness, the brightness of each pixel in the im- 



age plane is taken to be equal to the brightness of the corre- 
sponding point of the source plane to which it is mapped. Pixel 
mapped outside the circle of radius around the center of the 
source do not belong to the observed image. The observed mag- 
nification is then defined as the ratio between the total brightness 
of the image and that of the source. This gives us the flux magni- 
fication, which is equal to the size magnification because of the 
conservation of surface brightness. 

This procedure introduces an error related to the pixelization 
of the definition of the observed image. The relevance of this er- 
ror source was checked while running our simulations. In partic- 
ular we first performed the simulations for a smooth lens without 
substructures, and compared the magnification measured with 
the above procedure with the theoretical value of the magnifica- 
tion in the image centroid, given by the model. The results of 
this test are the following: 



^ = 0.018 ^ 

H iJ. 



0.027, 



(46) 



where Afi//j and cr(jj)/fi are the average relative offset and stan- 
dard deviation between the theoretical and measured magnifica- 
tions, respectively. The test was performed with a grid resolution 
such that the length of a pixel is a tenth of rg. The resulting typ- 
ical error in the magnification is of order of a few percent. This 
is much smaller than the ~ 30% error expected for a measure- 
ment of the magnification of an early-type galaxy with the use 
of the Fundamental Plane relation. Given this result, we adopted 
the same grid resolution for the actual simulations. 

In the simulations with substructures, a small fraction of the 
images display strong lensing features (multiple images, arcs, 
or rings). This is made possible by the fact that the deflection 
angle for a TSIS lens approaches a finite value as the distance 
to the center approaches zero, as can be seen from (l45i . For 
such strongly lensed images the observed magnification calcu- 
lated with the above method is very different from the magnifi- 
cation obtained with a smooth lens model. Such events are rather 
rare (typically a few per 1000 images). On the other hand, if such 
situation occurred in an actual observation it would be easily rec- 
ognized as a strong lensing feature. A highly distorted image, as 
an arc is, would not be suitable for a Fundamental Plane mea- 
surement. Therefore these images are removed from our analy- 
sis. 

At the same time we also took care not to make an excessive 
use of this procedure, for the following reason. The capability of 
a clump of a given mass to form arcs and multiple images de- 
pends on its internal structure. In our model we assumed TSIS 
mass profiles, but the real case is likely to be different from that. 
In fact, Meneghetti et al. (2003) showed that semi-analytic mod- 
els with NFW subhalos provide a better agreement with the ob- 
served statistics of arcs than models with SIS subclumps. Then, 
if in our simulation we observe an arc created by the presence of 
a massive subhalo it could be that the same clump with a more 
realistic internal structure would have not produced an arc but 
only a highly distorted image whose magnification could have 
been measured. For such images the value of the observed mag- 
nification is likely to differ significantly from the value inferred 
by assuming a smooth mass distribution. Thus by eliminating 
them from the analysis we would bias the results towards a better 
accordance between observed magnifications and smooth model 
magnifications. This problem is more relevant in the simulations 
in which clumps are treated as point masses, as they are more 
capable of producing arcs. 

On the other hand the most significant departures of the ob- 
served magnifications from the smooth case are for images in 
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the proximity of the most massive (m > IO^^Mq) subclumps. 
This is in agreement with the work by Meneghetti et al. (120071 
see discussion in 6.2.1). These are galaxy-scale objects and it is 
unlikely that they exist only in dark form. In other words, such 
substructures are likely to be easily identified by the presence of 
a luminous component. Therefore, for an image that lies in the 
proximity of one such object we should immediately suspect that 
part of the observed magnification is caused by the presence of 
this substructure and we would be warned of the bias caused by 
interpreting this data as tracer of the smooth component of the 
cluster Moreover, a number of gravitational lensing studies have 
been presented in which the lensing contribution of individual 
galaxies was incorporated into the analysis (see e.g. Natarajan et 
al. 120041) . Therefore, the same thing could be done for images 
of early-type galaxies close to cluster member galaxies, which is 
the situation discussed here. 

When multiple images are observed (such as the pair in Fig. 
|6ll, we take the magnification measurement obtained with the 
larger image only. In a realistic situation, if Fundamental Plane 
measurements are carried out on both the images, it would be rel- 
atively easy to label them as a pair and then to infer the presence 
of a substructure. We decided to adopt the more conservative 
approach in which the observer does not see the counterimage. 



1 " 

1 



I I 



Fig. 6. Multiple image system. This image (and others with sim- 
ilar properties) is not discarded, since the shape of the image 
on the left is not so strongly distorted and in general it may not 
be recognized as part of a multiple image system. In such cases 
only the larger image is considered for the analysis. 




Table 1. TSISs subhaloes 



Fig. 5. Strongly lensed image in the proximity of a massive sub- 
halo. Such images are removed from our analysis. This is done 
only after making sure that the arc feature would be recognizable 
under realistic observing conditions, i.e. that the arc-like shape 
would not be smeared out by a realistic PSF. In this example the 
larger arc is a few arcseconds long and could be easily identified. 



5.3. Results 

We studied image magnifications with three realizations of sim- 
ulated clusters, with mass fraction in substructure / of 0.05, 0. 10 
and 0.15. For each case we adopted two different models for the 
internal structures of the subclumps: TSIS and point mass. For 



/ 


Afi/lj 




N,-ej 


0.05 


0.005 


0.061 


4 


0.10 


0.014 


0.077 


6 


0.15 


0.018 


0.091 


20 



Notes. Diffefrences between magnifications in the presence and in the 
absence of substructure, for three cluster realizations. Clumps are mod- 
eled as TSISs. Each simulation run was performed by generating 10000 
images. A'„^ is the number of strongly lensed images that were rejected 
in each run. 



each case we generated A^^ = 10000 sources with the proce- 
dure described above. For each source, we studied the observed 
magnification with the presence of substructure, Hsuhs, and the 
magnification that would be observed with a smooth mass dis- 
tribution, Hsmooth- In Tables [1] and |2] we report the differences 
between the measured values of jisubs and Hsmooth- In particular, 
we introduce the average relative difference, defined as 



^ = — V 



1^ s lib s , / smooth , / 
f^smoolhj 



(47) 



and the expected mean relative error. 



1 ^ ^ I A'.?»/j.s',/ ^smooihj 
s ■ t \ l-^ smooth J 



(48) 



For completeness we also record the number N,-ej of rejected 
images in each simulation run. 

The quantity that is most relevant for our study is the disper- 
sion of observed magnification around the expected value, cr^. 
As expected, this quantity increases with increasing mass frac- 
tion in substructure. However, the values of this dispersion are 
somehow small. This is rather good news, because it means that 
the magnification of early-type galaxies is more sensitive to the 
smooth component of the mass distribution, which accounts for 
the bulk of the mass, and is therefore a quantity suited to con- 
strain the total mass of the lens. 
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Table 2. Point mass subhaloes 



/ 






A',., 


0.05 


-0.014 


0.052 


41 


0.10 


0.002 


0.063 


104 


0.15 


0.027 


0.075 


127 



Notes. Differences between magnifications in the presence and in the 
absence of substructure, for three cluster realizations. Clumps are mod- 
eled as point masses. Each simulation run was performed by generating 
10000 images. N„j is the number of strongly lensed images that were 
rejected in each run. 



not differ substantially, then it means that weak lensing errors do 
not play an important role and that the theoretical treatment of 
Sect.|4]can find applications in practical cases. 

Before facing the problem in full it is interesting to study 
how the errors in the weak lensing analysis alone influence the 
estimates of the total mass of the lens. In other words, we wish 
to clarify what is the typical error in the estimate of M in the 
hypothetical case of perfect Fundamental Plane measurements. 
This situation is simulated first and a realistic case in its full 
aspects is studied later 



The case of point mass substructures deserves further dis- 
cussion. As noted above, point masses are expected to pro- 
duce larger differences between the observed magnification and 
fJ^smoorh- From a first look at the results of Tables[r|and|2] it seems 
that the opposite situation is realized. However, in the simulation 
with point mass substructures the number of strongly lensed im- 
ages that is rejected is higher than in the TSIS case. This means 
that a significant fraction of the sources whose images are re- 
jected in the point mass case would be included in the analysis 
if lensed by a model with TSIS substructures. Presumably, the 
images of these sources are magnified by substructures and they 
contribute significantly to the value of the measured dispersion. 
This means that in the TSIS case the observed higher value of 
the dispersion is determined by a small number (less than 1 in 
100) of images, and by excluding them from the analysis we 
would obtain a dispersion not larger than the one observed in 
the point mass simulation. This result is in qualitative agreement 
with the work by Meneghetti et al. (I2007I I. They showed that 
in a more realistic cluster realization the probability of having 
a tangential-to-radial magnification (equivalent to observed axis 
ratio for circular sources) larger than 5 is less than 3% (see Fig. 7 
of the cited paper). They also claim that, in their case, substruc- 
tures account for 30% of the strong lensing cross section (the 
smoothed version of their cluster is still a critical cluster, unlike 
our case), meaning that the above percentage must be scaled ac- 
cordingly to be compared with our study. 

On the basis of these results we conclude that the magnifi- 
cation of early-type galaxies is little influenced by the presence 
of substructure. Substructure seems to play a significant role in 
the image formation only when present in large amounts, which 
is an unlikely scenario. Current estimates of the mass fraction 
in substructure based on numerical simulations give as typical 
values / < 0.10 (Tormen et al. 1998; Ghigna et al. 2000). These 
results give more significance to the technique of magnification 
measurement based on the use of the Fundamental Plane rela- 
tion, and set a sohd base for the adoption of the technique de- 
scribed in this paper for the purpose of solving the problem of 
the mass-sheet degeneracy. 



6. Testing the method 

The goal of the study presented in this section is to clarify to 
what extent the results obtained in Sect. |4] also hold in more re- 
alistic situations. 

We take a model cluster lens, apply the mass measurement 
technique to synthetic weak lensing and magnification data sim- 
ulated for this model, and then we analyze the results obtained. 
In particular, we compare the dispersion on the measurement of 
the mass obtained in these simulations with the dispersion ex- 
pected by considering the effects of Fundamental Plane mea- 
surements only, obtained from (l36T i: if the two quantities will 



6.1. Simulations 

In the following we will describe the ingredients necessary to set 
up the simulations. 



6.1 .1 . The lens model 

The model adopted to describe the lens is a Nonsingular 
Isothermal Sphere (NIS). The choice of a smooth model for 
the lens is suggested by the results of the analysis described in 
Sect.|5]on the effects of substructures on the magnification of the 
images of early-type galaxies. 

The parameters of the NIS model chosen for this test are the 
same adopted for the description of the smooth component of 
the lens model used in the simulations of Sect. |5] ?"c = 88 kpc 
and o",, = 1 156 km s"'. 



6.1.2. Weak lensing data 

For the simulation of weak lensing measurements we take in- 
spiration from the paper by Bradac et al. (2Q04j. The procedure 
adopted is the following. 

- We generate background galaxies with a uniform spatial dis- 
tribution in a 4 X 4 arcmin^ field of view (magnification ef- 
fects on the spatial distribution of images are neglected). 
Three different values of the number density are chosen: 
n = 30, 50, 70ai-cmin"l 

- Each galaxy is assigned a redshift, taken from the following 
distribution: 



(49) 



with zo = 2/3, as suggested by Brainerd et al. (|1996l l. This is 
a standard choice for weak lensing simulations. 
Each background galaxy is then assigned an intrinsic elUp- 
ticity drawn from a truncated Gaussian distribution: 



1 



27rcr2[l -exp(-l/2cr2)] 



exp{-|ef /2tr2), (50) 



with = 0.25. This is also a standard choice for weak 
lensing simulations (e.g. see Bartelmann et al. 119961 Seitz 
& Schneider [T997] l. 

For each galaxy, the resulting ellipticity is calculated from 
Eq. (H. 

An artificial measurement error is added, so that the resulting 
measured ellipticity e'" is 

e"' = £ + (51) 

where e^'''' is a random error generated from a Gaussian dis- 
tribution with dispersion cTg,.,. = 0.1. In adding the errors we 
ensured that |e'"| < 1. 
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The data are then processed with the finite-field inversion tech- 
nique of Seitz & Schneider (1 19971 1 . A grid in the image plane is 
defined. At each grid point j}, the average ellipticity (e) is es- 
timated from the observed ellipticities of background galaxies 
as 



(52) 



where Wiff) is a Gaussian with dispersion AO, such that nAO^ = 
12. 



6.1.3. Fundamental Plane measurements 

After the weak lensing reconstruction, which provides a model 
density map Ko(ff) up to the invariance transformation (IT2t . the 
simulation proceeds with the generation of A^fp = 20 early-type 
galaxies and the related Fundamental Plane measurements. This 
is done as follows. 

The position of each galaxy is generated randomly with a 
uniform distribution on the image plane (again, the efi'ects of 
magnification in the spatial probability distribution of images 
are not considered). Each galaxy is then assigned a redshift 
0.5 < z < 1.0. The upper limit reflects the redshift limit reached 
by the Fundamental Plane measurements carried out so far. The 
lower limit instead is set because the lensing signal for sources 
too close to the lens is too low. The simulation also requires a 
specification of the shape of the redshift distribution of the ob- 
served early-type galaxies, Pf^\z). This quantity depends on 
the intrinsic luminosity function of early-type galaxies, which 
in general varies with redshift, and also on the object selec- 
tion procedure. Given these uncertainties, for our simulations 
we adopted a uniform distribution to reduce the computational 
effort. It is also assumed that no error is introduced in the deter- 
mination of the individual redshifts. 



For each early-type galaxy, the quantity 



is then generated from a Gaussian distribution with 15% disper- 
sion, centered on the true value given by the model. This disper- 
sion should reflect the observed scatter in of the Fundamental 
Plane relation. In our case we adopted an optimistic estimate of 
this latter quantity. 

6.2. Results 

Under the observation conditions described above (A^pp = 20) , 
the expected dispersion on the estimate of the total mass of the 
lens , calculated from (l36l l and therefore ignoring weak lensing 
errors, is 



' exp 



o-(M) 
M 



0.21 



lo.l5/' 



(53) 



where crpp is the dispersion in of the Fundamental Plane, av- 
eraged over redshift. 

To better quantify the effects of weak lensing, the simula- 
tions have first been performed in the hypothetical case of per- 
fect magnification measurements (/c-^^^ = K{6j)), from which we 
obtained an estimate of the dispersion in the measurement of the 
total mass introduced by weak lensing errors only. The results 
of the simulations relative to this particular case are presented in 
16.2.11 while in l6.2.2l we report those obtained in a more realistic 
situation, in which the simulated magnification measurements 
have a 30% dispersion. 



Let us consider a mass density map kwl{0) obtained from 
weak lensing and transformed with (fT2T l to reproduce the exact 
value of the lens mass within the field of view. Since weak lens- 
ing reconstructions typically produce smoothed density maps, 
we expect Kwiiff) to underestimate the surface mass density in 
the central parts of the lens, and to be higher than the true value 
in the outer parts. This is indeed what is observed in the sim- 
ulations: in Fig. Qwe plot the difference kwl{0) - k(0) for an 
example case of weak lensing reconstruction. 

This property has some consequences on the process of fit- 
ting the density map to Fundamental Plane measurements. Since 
the expectation value of /c^^^* is approximately the true surface 
mass density, measurements close to the center will tend to give 
estimates k^^^^ higher than /(■wl(^) and therefore will bias the 
measurement of the total mass towards higher values (remem- 
ber that K-^i^iff) was defined as the density map that corresponds 
to the correct value of the mass). In contrast. Fundamental Plane 
measurements in parts of the lens plane where Ky^/i - k > Q will 
tend to bias the mass towards lower values. Thus, if we distribute 
the sources uniformly we expect that this potential source of bias 
can be overcome statistically. This effect is more evident when 
there is no error on k^^^^: in that case, the errors in the estimate of 
the total mass are only due to errors in the weak lensing analysis. 




Fig. 7. Contour plot of the difference Kwiiff) - k(0), where k(0) 
is the true surface mass density of the lens and kwl(0) is a weak 
lensing reconstruction of the lens yielding the same total mass. 
In the central region the reconstructed profile underestimates the 
surface mass density, while in a significant region of the image 
plane the opposite case occurs. 



6.2.1. Ideal case: perfect magnification measurements 

1000 simulation runs are performed for three different val- 
ues of the number density of background galaxies: n = 
30, 50, 70 arcmin"^, and two different numbers of Fundamental 
Plane measurements: A^pp = 1 and A^pp = 20. 
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Table 3. Weak lensing errors only. A^fp - 1 



n (arcmin ~) 


AM/M 


cr(M)/M 


30 


0.010 


0.32 


50 


0.002 


0.26 


70 


-0.007 


0.23 



Notes. Relative mean error and standard deviation of the measured 
value of the total mass of the lens, for three different values of num- 
ber density of background galaxies, n. 

Table 4. Weak lensing errors only. A^pp - 20 



n (arcmin ") 


AM/M 


cr(M)/M 


30 


0.010 


0.069 


50 


0.008 


0.060 


70 


0.006 


0.055 



In Tables [3] and |4] the relative mean eiTors AM/M = (M - 
M)/M and dispersion cr{M)/M obtained in these simulations are 
reported. As expected, the error in the determination of the total 
mass decreases as the number density of background galaxies is 
increased. However, the dependence of cr{M) on n is mild: this 
means that the number density of background galaxies is not a 
critical factor in the effectiveness of the method. It can also be 
seen that with a single magnification measurement the total mass 
is poorly constrained, as expected. This result shows that with a 
single local estimate of the surface mass density, no matter how 
accurate, it is difficult to break the mass-sheet degeneracy. 

On the other hand, in the case of A^pp - 20 Fundamental 
Plane measurements the dispersion is only a few percent. 

Another significant result is that there is practically no bias in 
the estimate of the total mass (AM/M ^ 0). This result is clearly 
a consequence of the assumption of a uniform spatial distribu- 
tion for the Fundamental Plane measurements. If we manage to 
pick a sufficient number of early-type galaxies more or less uni- 
formly distributed in the field of view there are good chances for 
the final measurement of the mass to be unbiased. This is a great 
advantage of the present technique with respect to the use of 
strong lensing information. Strong lensing features are typically 
limited to the central regions of clusters. Then, since the surface 
mass density in the central parts of a cluster obtained from weak 
lensing is typically underestimated, the inclusion of strong lens- 
ing data might lead to higher estimates of the total mass. With 
the present method this effect can be kept under control. 

In principle, substructures can introduce noise in the weak 
lensing signal as well. The nonparametric reconstruction method 
used here can recover features of the scale of the smooth- 
ing length or larger, as shown by Seitz & Schneider (11997) . 
However, King et al. ( 120011 ) and Clowe et al. ( 2004) showed that 
the effects of small scale substructures are of modest importance 
for weak lensing measurements. Therefore we can conclude that 
the results presented here do not depend on our choice of a sim- 
plified cluster model. 

Additional noise in the weak lensing measurements can be 
introduced by the presence of uncorrelated large scale structure 
along the line of sight. This issue was studied extensively by 
Hoekstra (20031 He found that the contribution of structures 
non associated with the cluster is important at large radii, but 
is negligible for the relatively small fields of view considered in 
our work (a few arcminutes). 



the lensing magnification by galaxy clusters 
Table 5. General case. A^pp - 20 



n (arcmin 


AM/M 


o-(M)/M 


30 


0.015 


0.22 


50 


0.003 


0.22 


70 


0.011 


0.21 



6.2.2. The realistic case 

Finally we report the results of the simulations performed in 
the most general case, in which realistic conditions for both 
the weak lensing and the Fundamental Plane measurements are 
simulated. In this case the simulations are performed for n = 
30, 50, 70 arcmin"^ and A^pp = 20. In Table|2]the results relative 
to 1000 runs are reported. 

The results confirm the picture emerged in 16.2.11 weak lens- 
ing errors play little role in determining the final dispersion on 
the estimate of the total mass. In fact, the observed dispersion is 
practically equal to the one obtained in ( l53l ) by assuming perfect 
weak lensing measurements. This test, although limited to a sin- 
gle lens model, confirms that it is possible, with a sufficient num- 
ber of Fundamental Plane measurements uniformly distributed 
on the image plane, to break the mass-sheet degeneracy, at least 
for Coma cluster-like lenses at intermediate redshift. 

Finally, to better illustrate the effect of a cluster on lensed 
early-type galaxies, we show in Fig.|8]a set of simulated FP mea- 
surements, compared with the Fundamental Plane relation ex- 
pected in the absence of lensing. A^pp - 20 objects in the redshift 
interval 0.5 < z < 1 .0 are generated and placed randomly behind 
the same cluster lens used for the previous simulations. In con- 
structing this plot, we assumed local values of the Fundamental 
Plane coefficients (a - 1.25, y6 = 0.32, y = -8.970 (J0rgensen 
et al. II996I recomputed for Ho - 65 km s ' Mpc ' by Treu et 
al. 120051 ). corrected for evolution following Treu et al. (120051 ): 
yiz) - 7(0) + 0.58z. The scatter of the FP is assumed to be 20% 
in Tg. The signature of the lensing signal can be clearly seen as 
an upward shift in the FP space. 



1.4 




aLogdo + (3 (SB)^ + 7 + 0.58^ 

Fig. 8. A^pp - 20 simulated FP measurements, plotted in FP 
space. The straight line is the Fundamental Plane expected in the 
absence of lensing, viewed edge-on. Positions in the FP space 
are coiTected for the evolution of the FP with redshift, following 
Treu et al. (I2005] ). 
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7. Conclusions and discussion 



In this paper we presented a new lensing-based method for the 
measurement of the mass of galaxy clusters. This method relies 
on the joint use of weak lensing data and magnification informa- 
tion, where the latter is obtained from Fundamental Plane mea- 
surements on background early-type galaxies. 

A statistical study of the method was carried out, and simu- 
lations were performed to test the importance of the presence of 
substructures and of errors in the weak lensing analysis for the 
success of the measurement. 

Our main conclusions are the following: 

- The most important quantity on which the effectiveness of 
the method depends is the mean surface mass density within 
the field of view of observation, k, while little role is played 
by the shape of the mass distribution. 

- Substructures contribute at most with a scatter of a few per- 
cent on individual magnification measurements. 

- Weak lensing errors introduce only a small dispersion on the 
final estimate of the total mass. 

On the basis of these results, we will now discuss which are 
the best lens candidates for an application of the present tech- 
nique. 

An important limit to the applicability of this method is the 
difficulty in performing Fundamental Plane measurements, since 
they require a significant amount of telescope time. A realistic 
number of Fundamental Plane measurements that can be per- 
formed in an observational campaign is ~ 20. Given this fact, 
we can fix A^fp = 20 and discuss which systems are best ana- 
lyzed with this number of magnification measurements. 

One of the most important factors in determining whether 
a cluster can be realistically studied with our technique or not 
is its redshift. The redshift must be be sufficiently high for the 
critical density to be low enough, to allow for higher values of 
K for a given physical surface mass density. On the other hand, 
the redshift must also be sufficiently low so that it is possible to 
find an acceptable number of early-type galaxies behind it for 
which Fundamental Plane measurements can be performed. The 
current observational capabilities and the lack of a calibration 
of the Fundamental Plane relation at very high redshifts set z ~ 
1 .0 as the highest redshift for which these measurements can be 
performed today. 

Bearing this in mind, we plot in Figure |9] the value of the 
critical density as a function of source redshift for three differ- 
ent values of the lens redshift. It can be seen that with a lens 
redshift Zd = 0.1 the resulting critical density is significantly 
higher than in the other cases at the source redshifts of interest 
and for this reason this case should be discarded. At the oppo- 
site end, for a lens redshift Zd - 0.4 the critical density is indeed 
the smallest for source redshifts higher than ~ 0.8, but the range 
of source redshifts for which the critical density is significantly 
small is limited to Zs > 0.6, and it shrinks rapidly for increasing 
Zd- On the basis of these simple considerations we conclude that 
a suitable redshift range for our lens cluster is Zd ~ 0.2 ^ 0.4. 
A thorough analysis of the problem would require a detailed 
knowledge of the redshift distribution of the observable early- 
type galaxies. 

Then we can ask which intrinsic physical characteristics a 
cluster should have to be efficiently probed with Fundamental 
Plane measurements. 

If we assume, conservatively, a mean dispersion of the 
Fundamental Plane relation of 20% in it can be shown from 
( |36] | that in order to obtain a mass estimate with a precision 
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Fig. 9. Critical density (in Kg m ^) as a function of source redshift for 
three different values of the lens redshift: Zd =0.1, 0.2, 0.4. 



of 20% or better with A^fp = 20 Fundamental Plane measure- 
ments with uniform distribution in redshift between z - 0.5 and 
z = 1 .0, the average surface mass density of the lens should be 
^>0.3. 

Clearly, the quantity k depends on the mass distribution of 
the lens but also on the size of the field of observation, which 
is set by the observer How critical is the choice of this latter 
quantity for determining a value k = 0.3? One can always re- 
strict the observations to the inner regions of a given cluster, to 
increase k. However, there is a limit set by the number density of 
background early-type galaxies for which Fundamental Plane 
measurements can be effectively performed. This number den- 
sity was estimated inJBL06jto be ~ 2 arcmin"^. Therefore, if we 
wish to find A^fp - 20 objects suitable for our purposes, we need 
to cover a field of view of 10 arcmin^. 

At this point we can study what is the minimum mass a clus- 
ter should have in order to have a mean surface mass density 
K - 0.3 within a circle of area A - lOarcmin^ (and radius 
6* - y/lO/n ^ l.Sarcmin). Fixing the lens redshift Zd = 0.3, 
this value of the surface mass density within the circle cor- 
responds to a value of the enclosed projected mass equal to 
Menc = 4.2 x 10^* Mq. For a NFW profile with concentration 
parameter c = 10, this value corresponds to a limiting virial mass 

14 



M200 = 8.0 X 10 



(54) 



and similar values hold for different values of the concentra- 
tion c. The quantity Mjqq is the minimum mass a NFW cluster 
should have in order to satisfy k > 0.3 within a circle of area 
A - lOarcmin^ in the sky. In that case the mass of the cluster 
within the circle can be measured with a 20% precision or better 
with 20 Fundamental Plane measurements and a weak lensing 
analysis. 

If we want to extend the analysis to the outer parts of our 
cluster, we must increase the number of Fundamental Plane mea- 
surements, to compensate for the reduced value of the average 
surface mass density k. 
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If we want to improve the sensitivity to less massive clus- 
ters, we must either increase the number of Fundamental Plane 
measurements for fixed aperture or choose a smaller aperture 
for fixed A^pp. In any of these cases, a higher number density 
of sources is needed, which cannot be obtained without an im- 
provement of the observational capabilities. This shows that the 
number density of observable sources is indeed a key factor and 
confirms that the lens redshift should not be too high. 

A possible alternative strategy is to select only the blighter 
objects, since it has been recognized that the FP scatter decreases 
with increasing mass (Treu et al. 2005 , van der Wei et al 2005 1: 
this would allow for a more precise estimate of the magnification 
for a given lens and fixed number of FP measurements. The fea- 
sibility of such an approach also depends critically on the num- 
ber density of background sources. 

On the basis of this discussion, we conclude that the value of 
Mjgy given by ( l54l i is an estimate of the minimum mass a cluster 
should have to allow for a mass measurement with the present 
method. The value is on the high side, but there are indeed 
many clusters that have observed values of M20() higher than this 
threshold. Well-studied examp les are A1689, A1703, A370, RX 
J1347-11 (see Broadhurst et al. l2008l for a review). All these sys- 
tems display strong lensing features that allow for a good esti- 
mate of the mass distribution in the inner (6 < 1 arcmin) regions 
of the clusters. One might think that the availability of strong 
lensing data would rule out the need for other observations in 
such clusters. Neveitheless, a great benefit would come from the 
addition of Fundamental Plane measurements as they could pro- 
vide important constraints on the mass distribution within a ra- 
dius at least 2 times larger 

In summary, the method presented in this paper is, because of 
its nonparametric form, a potentially powerful tool to break the 
mass-sheet degeneracy in lensing studies of clusters of galaxies. 
As shown in Sect. H) it allows also for a relatively easy estimate 
of the accuracy of the total mass measurement. 

In addition, the method can be extended by allowing for the 
inclusion of magnification measurements obtained from differ- 
ent means, for example from the observation of type la super- 
novae (see, e.g., Holz 2001; Goobar et al. 120021 Jonsson et al. 
1201 Oi l. The statistical framework developed in this paper can 
be applied with little effort to such more general situations, and 
therefore can be used as a reference framework to estimate the 
degree of precision of other methods that rely on the combina- 
tion of weak lensing and magnification measurements. 
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Appendix A: Probability distribution of the 
estimates of a^fp 

Here we will discuss the choice of a Gaussian probability distri- 
bution for the estimates a:fp of the surface mass density obtained 
by combining Fundamental Plane and weak lensing measure- 
ments. The observable quantity is the ratio r = rf^^ /r^f"'^ of the 
galaxy effective radius inferred from the Fundamental Plane re- 
lation to the observed (magnified) effective radius. Of these two 
quantities, the first, with its 15% dispersion, is by far the one 
with the larger uncertainty. Thus, if the probability distribution 
of the estimate of rf^^ is Gaussian, the ratio r will have a nearly 
Gaussian distribution as well. Let us assume this is the case. 

The surface mass density k^^^ at the image position can be 
expressed as a function of magnification /i, galaxy redshift z and 
average distortion (e), as in ( fTTl i. (fTST l, ( fT9] l. Let us assume that 
no error comes from the measurements of z and (e), in accor- 
dance with the hypotheses of Sect. |4] The only error source is 
then the measurement of the magnification. Note that, if /c*^^* 
depends linearly on the ratio r, then /c*^^* would also have a 
Gaussian probability distribution. However, the dependence of 
/f*-^' on magnification is more complicated than a simple linear 
relation. In fact, a:*^^' is given by 

^(FP) ^ -fc- y/b^-a{c-r'-) ^ ^^ ^^ 

a 

Now, if the condition 

- ac «i ar^ (A.2) 
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holds, then the square root can be approximated by 

VZ?^ - ac + ar^ ~ yfar (A.3) 

and (lA.lb becomes linear As a consequence, the Gaussianity of 
r translates into the Gaussianity of k''^^\ At this point, we need 
to demonstrate the validity of ( IA.2b . The ratio r is typically of 
order 1, unless in the case of large magnifications. Thus, it is 
sufficient to prove that - ac <s: a. We begin by noting that 
a ~ Zf[l - 0(\ (e) 1^)]. The quantity | (e) | is typically small, ex- 
cept in the proximity of critical curves, while the combination 
of weights that multiplies | (e) p in a is of order unity. Similarly, 
b ~ Z,[l - 0(1 <e) P)] and c ~ 1 - 0(1 <e> P). Then, the quan- 
tity b^ - ac is of order ZrO(\ (e) <K 1 . Since the right hand 
side of (IA.2I 1 is of order unity, (IA.2l i is satisfied. This implies that 
the dependence of a:'^^' on r is approximately linear Then, if the 
probability distribution of r is Gaussian, which is reasonable, the 
same applies to a:*^^'. 
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